DYNAMIC RESPONSE OF 
GFRP LAMINATES WITH EMBEDDED PZT 

PATCHES 


By 

Raja Yellamaraju 



DEPARTMENT OF AEROSPACE ENGINEERING 

Indian Institute of Technology Kanpur 


FEBRUARY, 2002 



DYNAMIC RESPONSE OF 
GFRP LAMINATES WITH EMBEDDED PZT 

PATCHES 


A Thesis Submitted 

In Partial Fulfillment of the Requirements 
for the Degree of 
Master of Technology 


by 

RAJA YELLAMARAJU 



to the 

DEPARTMENT OF AEROSPACE ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

FEBRUARY, 2002. 



: - 
??;c .:, 



137904 



CERTIFICATE 


It is certified that the work contained in the thesis entitled ‘^DYNAMIC 
RESPONSE OF GFRP LAMINATES WITH EMBEDDED PZT 
PATCHES", by Raja Yellamaraju, has been carried out under our su- 
pervision and that this work has not been submitted elsewhere for a degree. 


Prof S.Kamle 

Dept, of Aerospace Engineering 
I.I.T Kanpur 


y- * - "It 

Prof. N.G.R. Iyengar 

Dept, of Aerospace Engineering 

I.I.T. Kanpur 


FEBRUARY, 2002 



ACKNOWLED CEMENT 


I take this opportunity to express my immense gratitude and heartfelt thanks 
to Prof. Sudhir Kamle and Prof. N.G.R.Iyengar for their invaluable guid- 
ance, constant encouragement and advice at every stage of my thesis that 
made this work possible. 

I am extremely indebted to Mr.B. Prasad and Mr Lavendra Singh of Aero- 
structure lab for assisting me in fabricating various set ups required for this 
work. I would like to thank Mr.B.K Jain of ACMS department, Mr.Diwakar 
of material testing lab and Mr.M.N Mungole of materials science lab for al- 
lowing me to use various facilities for successfully conducting the experiments. 

My special regards to S.K.Muneshwar who helped me in experimental set 
up, Shanmukhi Gupta and to all my friends, who have made my stay at IIT 
Kanpur enjoyable. 

Last, but not the least, I must express my deepest sense of gratitude to my 
grandpa,parents,mentors(Mr M.S Raju and Mr Y.M.Raju) for their constant 
support and encouragement. 


Raju YRN. 



Abstract 


A smart laminated composite structure incorporates distributed actuators 
and sensors that are either surface bonded or embedded or both with the host 
structure and also includes onboard control and instrumentation. This work 
presents preparation of composite laminate from plain weave GFRP prepregs. 
The composite laminate was characterised by performing various tests to 
determine the four independent elastic constants. The smart glass/epoxy test 
specimen was manufactured with embedded PZT crystals. The test specimen 
was clamped on the fixture at one edge. The beam is connected to a magnetic 
shaker through a connecting rod at its other edge. The beam was subjected 
to a concentrated time varying(sinusoidal) load at its tip. The excitations 
of the beam were sensed by the embedded crystals. Accelerometer used as a 
sensor measures the response of the beam. The output of the crystals was 
amplified by using a charge amplifier and measured with the help of FFT 
Analyser. It was found that the response of the sensors were close to the 
theoretical estimates. Theoretical analysis was done by MSG Patran fern 
package. 
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Chapter 1 

INTRODUCTION 


1.1 GENERAL 

The word composite means, “consisting of two or more distinct phases”. Thus 
materials having two or more distinct constituents having distinct interface 
separating them are called composite materials. Composites can achieve 
certain physical properties not realizable by the constituent materials indi- 
vidually. These materials offer high specific strength and stiffness. Oriented 
fibrous composites also offer controlled anisotropy. The manufacture of com- 
posites requires comparatively low labor and generates less wastage in ad- 
dition to the ease of processing of complicated structural forms. Composite 
structures are destined to many future applications, especially in Aerospace 
and transport industry due to their properties such as good fatigue and cor- 
rosion resistance, low heat conductivity, good electrical insulation properties, 
cost effectiveness, etc. In the recent past, extensive research on the manufac- 
ture, characterization, fatigue and fracture of these materials have increased 
the confidence in the use of these materials and thus composites have replaced 
conventional materials in many applications. 

With progress in engineering design, increasing use is being made of 
lightweight high-strength materials. This is true irrespective of the field 
of their application whether it be buildings, bridges. Aerospace structures 
or mechanical structures. Spectacular progress has been witnessed in re- 
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cent years in a variety of scientific and engineering disciplines. In particu- 
lar, metallic and composite materials and structural technology, computer 
aided structural design and analysis, microelectronics, microprocessors, sen- 
sors, signal processing, etc have seen rapid advancement. The synergistic 
interaction of these fields has lead to the birth of the new trans disciplinary 
field of Smart/Intelligent Structures Technology. Some have viewed these 
developments as the “dawn of a new materials and structures age!”. ’Smart’, 
’Intelligent’ and ’Adaptive’ have all been used to describe and classify struc- 
tures which contain their own sensors, actuators and control capabilities. 
The “smartness” or “intelligence” can be defined at three levels, each of them 
defining a field of research, technology development and potential applica- 
tions: 


• A material or a structure is said to be “sensitive” when it includes 
sensors providing information concerning the material itself or its en- 
vironment. For instance the control of large and complex technological 
structures is a topic of rapidly growing interest. 

• A material will be “adaptable” if integrated actuators (or “active” ma- 
terials) can modify its characteristics. Such a material or structure will 
be adaptable only by the way of an externally determined action. 

• The combination of these two above mentioned properties results in 
an “adaptive” or “really smart” material which collects data related to 
the changes in its environment or in its own evolution or damage. It 
processes these collected data, and reacts through its ’actuators’ action. 

1.2 Sensors and Actuators 

Many types of sensors and actuators are being considered for introduction in 
smart systems. Typical sensors being used are of strain gauges, accelerome- 
ters, fiber optic sensors, piezoelectric films and piezoceramics. Sensors con- 
vert mechanical quantity (such as strain or displacement) into electric field. 
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Key factors for sensors are their sensitivity to strain or displacement , band- 
mdth and size. Other factors are temperature sensitivity, linearity, hystere- 
sis, electromagnetic compatibility, embeddability and associated electronics. 

Typical actuators are of piezoceramics, magnetostrictives, electrostric- 
tives and shape memory alloys. Piezoelectrics and electrostrictives are avail- 
able as ceramics, where as magnetostrictors and shape memory alloys are 
available as metal alloys. Most important performance parameters of actua- 
tors include maximum stroke or strain, maximum block force, stiffness and 
bandwidth. All these actuators directly convert actuating (electric/thermal) 
signals into actuation strain/displacement. 

There has been some significant work investigating shape control and 
damage control of structures but not nearly to the extent of active and pas- 
sive vibration control. Passive techniques utilize prior adjustment of the three 
design parameters - mass, stiffness and damping - to give a structure good 
dynamic response to an outside disturbance. These methods include tuning 
masses to place natural frequencies away from driving frequencies, vibration 
absorbers which create anti-resonance, and additional damping. Active con- 
trol techniques allow for a much larger range of operating conditions Present 
active control schemes require the use of applied external forces or motion at 
discrete locations to oppose the vibration and effectively dissipate the energy 
in the structure. A structure with the ability to change its stiffness or damp- 
ing in a controlled manner would offer an alternative to discrete external 
actuators. Smart structure applications are wide ranging from active shape 
control, vibration and noise control, improved damping, improved aeroelas- 
tic stability, to change of stress distribution. Many types of actuators and 
sensors are being considered, such as fiber optics, electro-rheological fluids, 
piezoelectric materials, shape memory alloys, electrostrictive materials, mag- 
netostrictive materials, etc. These materials can be embedded in composite 
laminates and can be used as integral sensors or actuators. 
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1.3 Smart materials 

The design of smart structure is significantly dependent on the appropriate 
selection of smart materials. The available smart materials can be classified 
into five major groups based on the characteristic working principles[l-3]. 
They are piezoelectric, magnetostrictive, shape memory, electro-rheological 
and magneto-rheological materials. Piezoelectric materials develop charge if 
deformed by the mechanical stress and by a converse effect, it deforms due 
to the application of electric field. A similar direct and converse relationship 
between mechanical and magnetic field exists for magnetostrictive material. 
Shape memory alloys have special features of memorizing a shape stretching, 
bending or twisting as a function of temperature and recovering that shape at 
different temperatures. They deform due to a phase change from martensite 
to austenite state. The phase transition is dependent on both stress and 
temperature such that by changing any one of them, volume change can 
be initiated. Electro-rheological fluids are a class of specially formulated 
colloidal suspensions, which undergo a change in the resistance to flow due 
to an electric field. Magneto-rheological fluids can be very effectively used 
for vibration suppression by building actuators out of this material. 

The current generation of smart structures featuring piezoelectric ma- 
terials is generally synthesized with polymeric fibrous composite laminates 
which readily accommodate embedded piezoelectric actuators and sensors. 
Any external force, applied on the beam will set vibrations and cause defor- 
mations in the beam. These deformations will cause stresses and strains in 
the beam. If the beam is embedded with smart crystals, these vibrations can 
be dampened quickly and deformations and resulting effects can be controlled 
to greater extent. 

Among the various materials, the piezoelectric materials have proved to 
be the best for sensors and actuators in smart structures. Piezoelectrics un- 
dergo surface elongation when an electric field is applied across them and 
produce voltage when surface strain is applied and thus can be used both 
as actuators and sensors. Various commonly used actuation materials are 
the piezoelectric materials - Lead Zirconate Titanate (PZT), and Polyvinyli- 
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denefluoride (PVDF); ceramic electrostrictive material - Lead Magnesium 
Niobate (PMN); Terfenol, a rare earth magnetic like material and Nitinol 
(nickel titanium alloy) available in the form of wires or sheets. These sub- 
stances have the ability to change the shape, natural frequency, damping 
and other mechanical characteristics of intelligent materials in response to 
a change in environment. When integrated into a structure either through 
embedding or surface bonding , they apply localized strains through which 
the deformation of the structural elements could be controlled. Due to their 
small size, a large number of piezoelectric actuator elements could be used 
without greatly increasing the mass of the structure. 
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LITERATURE REVIEW 


The problem of exerting control on large or precision flexible structure has, 
in recent years, been the center of attraction of an exhaustive research ef- 
fort through out the Aerospace engineering community [4]. Many presently 
planned and future missions such as structural health monitoring, vibration 
control, airfoil shape control, noise control, control surface actuation, smart 
radar, etc. need for pointing accuracy for better performance execution. 
World wide, many researchers have paid their attention towards the devel- 
opment of intelligent structures and have shown encouraging results in this 
field. The smart structure can help to reduce amplitude of vibration, bring 
down the overall mass, reduce drag on the wing and noise inside a fuselage, 
to name a few. They can be used on all Aerospace vehicles [5]. 

A detailed literature survey of the issues related to properties of PZT 
ceramics and their use as actuators for vibration control is not attempted in 
this chapter because of the extensive activity in this field. However, typical 
references are mentioned which provide some measure of background for the 
various issues of concern. 

The actuators in a smart structure provide the mechanism for the struc- 
ture to adapt to its surrounding by suppressing vibrations or changing the 
structure’s shape. This requires an actuation authority over a broad band- 
width spectrum and over a wide range of displacement amplitudes (from 
tenths of micrometers to milli/centi meters). 

Hagood et al [6] have shown that piezoelectric crystals can readily be em- 
bedded in composite structures without much affecting the global stiffness of 
the substructure and increasing the mass of the total system, to maximize the 
energy dissipation in the structure and exhibiting control on the vibrations. 
Another feature of their work is the development of suitable technique of 
embedding piezoelectric crystals in a composite laminate. They have shown 
that piezoelectrics are uniquely suited as elements of highly distributed sen- 
sors and actuators in smart structures. 

Crawley and Luis [7] have investigated the use of piezoelectric actuators 
to excite steady-sate resonant vibrations in simple structural members such 
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as catilevered beams. The response of the specimens were found in good 
agreement with those predicted by the analytical models. 

The feasibility of adopting piezoelectric films to detect damages in com- 
posites has been discussed by Shen et al [8]. In order to study the possi- 
bility of detecting and monitoring damage in composite laminates through 
piezoelectric film sensors, three different laminates in the presence of various 
damage were investigated and influence of damage on damage monitoring 
sensors was specifically analyzed. 

Zhou et al [9] have presented a dynamic analytical approach for the de- 
sign and integration of active piezoceramic (PZT) patch elements locally 
coupled with host structures. Several critical design issues have been ad- 
dressed like determination of the actuator dynamic outputs, the prediction 
of energy conversion efficiency, the estimation of system power requirement 
and the limitation of induced alternate peak stress. Both the mechanical 
stress behaviour and the thermal stress characteristics of the PZT patch ele- 
ments were investigated. The attention in parametric design was directed to 
the thickness and location of the elements. 

Batra et al [10] have shown the determination of the optimum location 
of a given rectangular piezoelectric actuator that will require the minimum 
voltage to annul the deflection of a simply supported rectangular elastic plate 
vibrating near one of its fundamental frequencies. 

Lee [11] presented theory of laminated piezoelectric plates for the design of 
distributed sensors/ actuators. In his theory, the piezoelectric phenomenon to 
effect distributed control and sensing of bending, torsion, shearing, shrinking 
and stretching of a flexible plate has been developed. The theory is capable of 
modeling the electromechanical (actuating) and mechanoelectrical (sensing) 
behaviour of a laminate. 

Padma Akella et al [12] studied the modeling and control issues related 
to smart structures bonded with piezoelectric sensors and actuators. They 
applied Hamilton’s principle to obtain a linearized equation of motion. The 
natural modes are then found by solving an eigenvalue problem. Herman 
Shen [13] developed a one - dimensional theory for modeling the analysis 
of beams containing piezoelectric sensors and actuators. The equations of 
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motion and associated boundary conditions are derived for the vibrations of 
piezoelectrically sensored/actuated beams. 

Chandrashekhara and Agarwal [14] presented a finite element formulation 
for modeling the behaviour of laminated composites with integrated piezo- 
electric sensors and actuators. This model is valid for both continuous and 
segmented piezoelectric elements that can be either surface bonded or em- 
bedded in the laminated plate. The formulation is based on the first order 
shear deformation theory, which is applicable for both thin and moderately 
thick plates. Liu et al [15] presented a finite element formulation to model the 
dynamic as well as static response of laminated composite plates containing 
integrated piezoelectric sensors and actuators subjected to both mechanical 
and electrical loadings. 


8 



Present Work 


Investigations on dynamics of smart structures have been drawing attention 
of many researchers worldwide. Smart materials, may be defined as mate- 
rials that posses the adaptive capabilities to external stimuli such as load 
or environment with inherent intelligence. The intelligence of the material 
could perhaps be programmed by material composition, processing, defect 
and microstructure or conditioning to adapt in a controlled manner to vari- 
ous levels of stimulus. Piezoelectric ceramic materials are typically employed 
as sensors and actuators. The objective of the current work are to perform 
various tests on the composite laminate to get the four independent constants 
and to study the vibration response of composite beam embedded with PZT 
crv'stals. 

The various tests employed for characterising the laminate are described 
in Chapter 2. Formulation for free vibration response of the plate is given 
in Chapter 3. Theoretical analysis of beam vibration response is described 
in Chapter 4. Chapter 5 deals with the results, discussions, conclusions and 
suggestions for future work. 
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Chapter 2 

PREPARATION OF LAMINATE 

2.1 Basic Raw Materials: 

The basic raw materials for the preparation of the laminate are the reinforcing 
material (glass woven fabric) and the matrix material (epoxy mixture). 

2.1.1 Glass woven fabric 

Woven fabrics are one of the most widely used reinforcing materials. Fab- 
rics typically offer flexibility in fabrication techniques, but at a higher cost. 
Woven reinforcing fabrics are made by interlacing individual filaments, ends 
(untwisted fiber bundles), yarns (twisted fiber bundles), and rovings. Fabric 
composites have mechanical properties similar to those of laminates made 
from orthogonal unidirectional layers. Fabrics can be woven into many dif- 
ferent types of weave patterns, widths and thicknesses. A variety of weave 
patterns can be used to interlace the warp and weft (filling) yarns to form 
a stable fabric. The major fabric styles are plain weaves, twills, satins and 
woven rovings. The plain weave, which interfaces one warp yarn over and 
under the weft (or filling) yarn, demonstrates the greatest degree of stability 
with respect to yarn slippage and fabric distortion (stability). Yarn count 
and content, however, also contribute to fabric stability. In the present work, 
plain weave fabric is used as the reinforcing material. 
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2.1.2 Epoxy Resin System 

One of the most versatile properties of epoxy resins is their ability to trans- 
form readily from liquid state to tough, hard thermo-set solids. This harden- 
ing is accomplished by the addition of a chemically active reagent known as 
curing agent or hardener (sometimes requires heat and pressure). Accelerom- 
eters are sometimes added in catalytic amount to reduce the cure time. The 
specific resin-curing system will, to a very considerable extent, determine the 
ultimate properties of the laminate. 

2.2 Composite Laminate for Dynamic Test 

2.2.1 Specimen Preparation 

Woven glass fiber cloth was used for making specimens for Dynamic and 
Tensile tests. The cloth was cut slightly larger than the required dimension 
into 19 layers. These are weighed to calculate the amount of epoxy required. 
10% extra amount of the epoxy was added to make up the loss during lay 
up. Hand lay up technique was used to prepare the specimen. After lay up, 
the specimen was covered on top with a mylar sheet and a flat glass plate 
is placed on it. Weights were put on the glass surface and the extra matrix 
applied during the lay up was removed. Care should be taken while placing 
the weights, so that no slippage occurs between the layers. It was cured 
at room temperature for about 24 hours. The specimen was removed after 
curing and was cured in the oven at 100-120° C for 2 to 3 hours. Final shape 
of the specimen was obtained by discarding the ends of the cured specimens. 
The size of the test specimen was 260 mm by 80 mm by 3.40 mm. 

2.2.2 Embedding PZT Crystals into the Test Specimen 

To accommodate the PZT disc’s in the laminate, fibers from that particular 
location should be removed. The size of the crystal was 10 mm by 10 mm by 
0.25 mm. Lead wire of diameter 0.2 mm was soldered on both faces of the 
crystal. The holes are made 1 mm greater than the size of the PZT crystals. 
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Two crystals are used and each is placed at equi-distance from the mid layer 
of the laminate at a distance of one - fourth the length from one edge and 
at the centre of plate width. The specimen with the embedded crystals is 
shown in the Fig.2.1. 

2.3 Characterization Of Composite Laminate 

2.3.1 Tensile Test 

Tensile modulus of the composite laminate was first determined by the ten- 
sile test, performed on the specimens prepared as per the ASTM D3039-76 
[16] standard specifications. The properties deduced from the tension tests 
on composite materials are effective (average) properties. The test method 
(ASTM D3039-76) applies to unidirectional laminates but can also be per- 
formed on multi-directional laminates, woven fabrics, or discontinuous fiber 
composites. 

Specimen Preparation 

The most commonly used specimen geometries are the dog-bone specimen 
and flat specimen with end tabs. In the present work, flat specimens are used. 
The standard specimen used for the determination of the Young’s Modulus 
in the longitudinal direction of the composite laminate is shown in figure. 
Five specimens of same dimensions are used for the test. 

Woven glass cloth was cut into number of pieces required for the given 
thickness. The cloth was cut into the rectangular pieces of the size 230 mm 
by 150 mm. 18 plies were used to achieve the required thickness. Resin is 
heated to 70 - 80 °C and cooled to room temperature. Once it is cooled, 9-10 
parts by weight of the hardener is mixed thoroughly until even distribution 
was achieved. Once it is mixed, the prepared mixture is used without delay 
since the pot life of the mixture is about 0.5 - 1 hour. The Resin - Hardener 
ratio is given in Table 2.1 

Two flat glass plates of size 500 mm by 500 mm by 5 mm are used for 
laying the specimen. One plate is placed on a smooth surface with a mylar 
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sheet of same area on it. Hand lay up technique is used for the specimen 
preparation. The epoxy mixture is applied on each layer of the lamina with 
a brush and stacked evenly. Once all the plies are stacked, one more mylar 
sheet is placed on it. The other glass is placed on the top this mylar sheet. 
Weights are placed on the top surface of the glass for pressure. Extra resin 
applied during lamination could be drained out when pressure is applied over 
the laminate. Curing was carried out at room temperature for about 24 hours 
and then semi-cured laminate was transferred to the oven, where it was held 
at 100 - 110 °C for about 1-2 hours. Finally, the end parts are discarded 
and the five specimens of required size are cut from the laminate. The tensile 
test specimen is shown in Fig. 2. 2. 

Points to be cared during specimen preparations are: 

1. All the surfaces must be flat, any irregularity should be avoided. 

2. Resin, if needed, should be heated to 70 - 80 "C to avoid moisture 
trapping before mixing with hardener. 

3. PVA solution must be applied over the surfaces of every part. 

4. Pressure should be applied uniformly on the laminate to reduce the 
non-uniformity in thickness. 

5. Specimen must be symmetric with respect to centre line. 

Testing Procedure 

The tensile test of all the five specimens was carried out on a MTS - 810 
machine. The tests were done in the load control mode. The specimen was 
clamped between the hydraulic grips and aligned correctly. Load was applied 
at the rate of 0.4 ton/min. An extensometer (MTS 632 - 25C - 20) was used 
to measure the strain. 

A typical load vs strain curve for the specimen is shown in Fig.2.3. The 
data obtained from the test was analysed to estimate the strength and elastic 
modulus of the laminate. 
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2.3.2 Burn out test 

Burn out test was carried out to calculate the exact fiber fraction in the 
laminate and finding out the various properties of the laminate by applying 
the rule of mixtures. The significance of this method [17] is it can be used 
to obtain the ignition loss of a cured reinforced resin sample. If only glass 
fabric is used as the reinforcement of an organic resin that is completely 
decomposed to volatile materials under the conditions of this test and the 
small amount of volatiles (water, residue solvent) that may be present is 
ignored, the ignition loss can be considered to be the resin content of the 
samples. This method does not provide a measure for resin content of the 
samples containing reinforce materials that lose weight under the conditions 
of the test or containing resin that does not decompose to volatiles materials 
released by ignition. 

Apparatus 

Crucible ceramic crucibles were used. Electric furnace, variac controlled elec- 
tric furnace capable of maintaining 1350 °C was used for this test. 

Test Specimen 

For the sake of homogeneity, samples were cut from the specimens that had 
been used for mechanical properties determination. The fractured areas were 
removed and three specimens of the sizes given in Table 2.2 were prepared. 
Thickness of all the specimens was 3.10 mm. 

Procedure 

The crucibles were heated to 500 to 600° C for 15 minutes or more, cooled 
to room temperature in a desiccator and weighed to the nearest 0.1 mg. 
The composite specimens were also weighed to the nearest 0.1 mg. Then 
they were placed in the respective crucibles. Weight of all the crucibles, the 
specimens and the weight of the crucible with specimen before and after test 
axe summarized in the Table 2.3 
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Crucibles and the contents were placed in the furnace and the temperature 
was raised up to 600° C. they were held at this temperature for 6 hours, the 
time until all carbonaceous material has disappeared. They were then cooled 
to room temperature in the desiccator and weighed to the nearest 0.1 mg. 
Since weight of the fiber and matrix and density of both the fiber and matrix 
are known, one can easily calculate the volume fraction of fiber. The volume 
fraction of the fiber was found to be 0.46 on average. The difference in fiber 
fraction of the various samples was within 0.8%. The density of the laminate 
was 1 643 gm/cm3. 

2.3.3 Inplane Shear Test 

Test Method 

Inplane shear properties of a composite can be determined by conducting a 
tension test on a [±45°]^ glass/epoxy fabric laminate The test results are 
interpreted by carrying out stress analysis of the laminate. For an applied 
stress Gxi we have, 


O'x 

nT = Y 


and 


F 


cr~ — 


(bxd) 

The inplane shear strength is, therefore. 


t’lt — 


2 X (b X d) 


where F is the failure load in Kilo Newton (KN), b and d are the cross- 
sectional dimensions of the specimen in mm 

The lamina strains are related to the laminate strains as follows: 


7Lr 


= - ^y) 
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The inplane shear modulus is, therefore, 


(jtlT = 

iLT 

F 

2 X (6 X d) X (ej - e") 

The Mohr circle for stresses and strains are shown in Fig.2.4 and Fig.2.5 
respectively. 

Specimen dimensions 

The geometry of the [±45°]4 j tensile off-axis is shown in Fig.2.6. Two spec- 
imens were used for testing. The size of the specimens were (190 x 25 x 
2.65)mm. The end tabs were of aluminium. The size of end tabs were 
(40 X 25 X 2) mm. 

Test Procedure 

For characterizing the material, tension tests were performed. The test was 
carried out on Table Top Tensile Testing Machine (Hounsfield Test Machine 
- W series: Maximum load - 20 KN) . The specimen is first aligned in the 
grips and tightened in place. The specimen was then loaded at a constant 
rate of 1.5 mm/min. The applied load was measured from the display menu 
on the machine. In order to determine the specimen strains, strain gauges 
were mounted on the specimen and monitored during the test. Strain gauges 
with resistance of 120 ±0.2 Q, were used. Two gauges were applied to each 
specimen in the longitudinal and transverse to the loading direction. The 
gauges were located close to the geometric centre of the specimen. The gauge 
reading was measured with the help of a multi channel strain indicator. A 
stress-strain diagram is drawn from the test data and shown in Fig.2.7. From 
the data, the value of inplane shear modulus was calculated. 
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Table 2.1: RESIN HARDNER RATIO 



Parts by weight 

Araldite CY 230 

100 

Hardner HY 951 

9- 10 


Table 2.2: BURN OUT TEST SPECIMENS SIZE 


S.No 

Size of the Specimen (mm x mm) 

1 

2.00 X 2.15 

2 

2.05 X 2.15 

3 

2.00 X 2.20 


Table 2.3: WEIGHT OF THE SPECIMEN BEFORE AND AFTER BURN 


OUT TEST 


S.No 

Specimens 
wt. (gm) 

Crucibles 
wt. (gm) 

Crucibles wt. with 
specimen (gm) 
(before test) 

Crucibles wt. with 
specimen (gm) 
(after test) 

1 

2.1459 

45.9964 

48.1423 

47.2033 

2 

2.2437 

37.3942 

39.6379 

38.6390 

3 

2.2885 

41.8910 

44.1795 

43.1620 


2.3.4 Determination of Poisson’s Ratio 

The size of the specimen used was 160 mm by 25 mm by 3.10 mm. Strain 
gauges were applied in longitudinal and transverse direction at equi-distant 
from centre of the specimen. The specimen is aligned in place between the 
grips of the machine. The strain gauge readings were recorded by a multi- 
channel strain indicator. Poisson ratio was calculated by dividing lateral 
strain by longitudinal strain. 
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Aluminium End Tab 



Tensile Test Specimen 

^ All Dimensions are m mm 

Figure 2.2: TENSILE TEST SPECIMEN 
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Gauge Length = 50 mm 

Load 100% ( 10 V= 10,000 kg) 

Strain 50% (10 V = 12.5 mm displacement or 0.250 strain) 

Stroke 10% (10 V = 10 mm) 


Figure 2.3: A TYPICAL LOAD VS STRAIN CURVE FOR TENSILE TEST 
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Figure 2.4; MOHR CIRCLE DIAGRAM FOR STRESSES IN INPLANE 
SHEAR TEST 
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Figure 2.5: MOHR CIRCLE DIAGRAM FOR STRAINS IN INPLANE 
SHEAR TEST 
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40 

AlumimumEndTab 



Shear Test Specimen AlDimensions are in mm 


Figure 2.6: [± 45 °] TEST SPECIMEN FOR INPLANE SHEAR TEST 
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Figure 2.7: A TYPICAL LOAD VS STRAIN CURVE FOR INPLANE 
SHEAR TEST 
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Chapter 3 

FORMULATION FOR FREE 
VIBRATION 


The problem is formulated for a fiber - reinforced laminated composite rect- 
^ angular plate of length a, width b, and thickness h composed of orthotropic 

■ ’ii layers of thickness hi with fibers oriented at angles 0 and -9 as shown in figure 

i 3.1. 

Certain assumptions made for the finite element formulation were: 

' 1. The area and mass of the actuator embedded in the plate is negligible 

‘ ■ compared to that of the plate. 

“ . , 2. Equations for actuator and sensor were not considered. 

3.1 Strain-Displacement Relations 

^ ' ’XTThe displacement field based on a first - order shear deformation theory is 

( r. rqriven by (Yang et al. [18]) 

u(x, y, z, t) = u°{x, y, t) + z'lp^ix, y, t) 


v{x, y, z, t) = v°{x, y, t) + zipy{x, y, t) (3.1) 
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w{x,y,z,t) = w°{x,y,t) 

where u, v, and w are the displacement components in the x, y, and z 
directions, respectively, t is the time, u°, v°, and are the in - plane and 
transverse displacements of a point (x,y) on the mid-plane, respectively, V’c 
and tpy are the rotations of the normal to the mid-plane about the y and x 
axes, respectively. 

Equation 3.1 can be expressed in matrix notation as 

{«} = [G]{u} (3.2) 

where 

{m} = (u, (3.3) 


{u} = 


(3.4) 


[G] = 


10 0^0 
0 1 0 0 z 

0 0 10 0 


(3.5) 


Substituting Equation 3.1 into infinitesimal strain relations, we can write 


{£} = {£“} + Z{X} 


(3.6) 


where 


{^} — (^i) ^z) T'xj/) Txzi Tyz) 


(3.7) 


0 , 7 ^, 7 °^, 
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' du° dv° ^ du° dv° dw° , dw° . '' 
■> O5 + “> + '^x, ^7- + V'y 


dx' dy' ' dy dx ’ dx 


dy 


(3.8) 


[x] — (Xz) Xy) 0) Xzyj 0) 0) 


' dipx 

dx 


d^y d^y ^ 

’ ’ -- dx ’ ’ > 


(3.9) 


dy dy 

The generalized strain vector {e} corresponding to the mid-plane can be 
defined as 


{e} = (el, e°, jly, 7 "^, 7 °^, Xi, Xy, Xzy)^ (3.10) 

which can be written as 


{6} = [L] {«} 


where [L] is a derivative operator matrix defined by 


where 


[L] = 


40000 
04000 
440 0 0 
00410 
00401 
00040 
00004 

0 0 0 dy dx 



and 


4 — 


1 . 

dy 


(3.11) 


(3.12) 


(3.13) 
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3.2 Stress strain relations for an orthotropic 
lamina 


Stress strain relations for a lamina in the material coordinate 
axis, whose fibers are oriented at an angle a with reference to the x — axis 
is given as [19-21] 


f N 

0-1 


Qn 

<5i2 

0 

0 

0 


’’ ^ 

ex 

0-2 


Q 21 

Q 22 

0 

0 

0 


€2 

ri2 

> = 

0 

0 

Qee 

0 

0 

< 

7l2 ' 

TlZ 


0 

0 

0 

Q44 

0 


7X3 

T 23 


0 

0 

0 

0 

Q 55 


723 


(3.14) 


where (Ji, 0-2, ci, andea are stresses and strains in the direction parallel 
and perpendicular to the fiber direction,ri2, tis, T23, 712, 713, andT2z are shear 
stresses and shear strains in the respective directions. 

The stiffness coefficients Qn , Q22, , e.tc., are defined in terms of material 

properties as 


Qn = 


El 

(1 — ^'I2i^2l)’ 


Qi2 — 


(1 - 1/121^21)' 


Q2I — Qi2 , 



E2 

(1 — Vi2l'2iy 


(3.15) 


Qm — Gi2i 
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<344 — G^13) 


Q55 — ^23- 


3.3 Laminate Constitutive Relations 


The stresses and strains in x, y and 2 directions are obtained 
by transforming the relation given in Equation 3.14. The transformed stress- 
strain relations for a lamina are given as 




Qii 

Qi2 

Qi6 

0 

0 


( 

^x 

(Jy 


Q 2 \ 

Q22 

Q26 

0 

0 


Cy 

^ '^xy 

> = 

Q\& 

Q26 

Qee 

0 

0 

< 

Ixy ^ 

'^xz 


0 

0 

0 

Q44 

Q45 


^xz 


k 

0 

0 

0 

Q54 

Q55 _ 

k 



where 

Qn = Qnrn^ + 2(<3i2 + 2Qee)m'^n^ + <322?^^, 


(3.16) 


<3i 2 = (<3n + Q 22 - 4(366)m^n^ -I- Quim^ + n^), 


Q 22 " Qn'^^ + 2((3i2 + 2 QQ^)nT?TT^ -h <322^^? 


Q\& — (^11 ~ ^12 — 2 Qm)m^n + {Q12 — Q22 + 2 Q^^)mn^ ^ 


Q 26 — (Qii ~ Qi2 “ 2QQQ)mn^ -f ((3i2 “ Q 22 + 2<366)^^^) 
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^66 — (Qii + Q22 — 2Qi2 — 2 Q^s)‘^'^ + 

<3 44 = Q 44771^ + Qss^^) 

Q45 = iQu - Q55)mn, 

Q55 ~ Qss"^ “(■ ^44^ 

where 

m = cos CK and n = sin a 

Using the above lamina constitutive equations, the laminate constitutive 
relations can be written as 

{W} = [5] {e} (3.17) 

where {jv| is the stress and moment resultant vector corresponding to 
the mid-plane defined by 

{iv} = (A^^, Ny, N^y, My, M ^y , Q^, Q y^ (3.18) 

rhl2 

{Nx^ ^xyi Qxt Qy) ^ I {f^xi^y'}'Txyi'^xz’t'^yz)dz 

J ’—hl2 


rh/2 

{Mx, My, Mxy) — I (<^X) ^yi '^xy) zdz 
J —h/2 
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^11 

4-12 

4i6 

Bii 

Bi2 

Bie 

0 

0 

■421 

422 

426 

B21 

B22 

B26 

0 

0 

Aiq 

426 

466 

Bie 

B26 

Bee 

0 

0 

Bn 

Bi2 

Bi6 

Dn 

Di 2 

Die 

0 

0 

B21 

B22 

B26 

D21 

D22 

D 2 e 

0 

0 

Bi6 

B26 

Bee 

Die 

-D26 

Dee 

0 

0 

0 

0 

0 

0 

0 

0 

444 

445 

0 

0 

0 

0 

0 

0 

454 

455 


By, By) = / Qy (1, 2, {h j = 1, 2, 6) 

fc=l 

= E / ' = 4,5) 

fc=i 


(3.19) 


(3.20) 


3.4 Energy equations 


Hamilton’s variational principle is used to derive the laminate 
equations of motion. The mathematical statement of Hamilton’s principle in 
the absence of damping can be written as 


r (5K -5U + 6W) dt = 0 
Jti 


(3.21) 


where K is the kinetic energy, U is the strain energy, and W is the 
potential energy of prescribed surface tractions. 

The kinetic energy term is 


5K= [ p{5uf {u]dV (3.22) 

J 

where 

{tt} = d {m} /dt is the velocity vector of any generic point (x, y, z) in 
space, p is the mass density of the material, and V is the volume of the 
plate. 
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Differentiating Equation 3.2 with respect to time, substituting into Equa- 
tion 3.22, and carrying out the integration through the thickness of the lam- 
inate, we get 

6K = J^{5 [m] {^} dA (3.23) 

where 

[»] = 

/i 0 0 /2 0 ' 

0 /i 0 0 /2 

0 0 /i 0 0 (3.24) 

/2 0 0 J3 0 

0 /2 0 0 /3 , 

where Ii, I 2 , and Iz are the normal, coupled normal - rotary and rotary inertia 
coefficients defined by 

{h, h, h) = Y^ [ Pk ( 1 , 2 ;, z'^) dz 

fc=i 

The strain energy term is 

6U= [ {5ef{a}dV (3.25) 

J \A 

Integrating above equation through the laminate thickness, gives 

5U = {5lf {]V} dA (3.26) 

Substituting Equation 3.17 into the above equation, gives 

6U= [ {(5e}^ [d] {e}dA (3.27) 

jA. 

For free vibration analysis, the potential energy 
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6W = 0 (3.28) 

Substituting Equations 3.23, 3.27, and 3.28 into Equation 3.21, we get 




dAdt = 0 


( 3 . 29 ) 


^t2 

Itl JA L' 

Integrating the first term by parts, the above equation can be written as 




3.5 Finite Element Formulation 

The finite element formulation is based on the first order shear deforma- 
tion theory. The plate is divided into a number of isoparametric four- 
noded quadrilateral elements with five degrees of freedom(DOF) per node 
{u°, v°, w°, V’l, tpy) ■ Linear interpolation functions for all variables over each 
element was taken. Accordingly, the displacement vector at any node within 
the element is expressed as 


{lf(a;,2/,t)} = ^[A;®(a;,y)]{u^(t)} (3.31) 

1=1 

where NN is the number of nodes per element and 

ra = -w.' (/] 

{Hf} = Vra.V’j.)’’ 

Here, are the element shape functions, [J] is a 5 x 5 identity ma- 
trix, and u", v°, w°, ■tjjxi and are the nodal values of u°, v°, w°, tpx and ipy, 
respectively . The superscript e denotes the parameter at element level. 
Equation 3.31 can also be written as 
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{tf} = [^f'] {a'} 


(3,32) 


where 


[iV) = [(Aff 1 . [W|| , 




{a®} — 7 7 

Substituting Equation 3.32 into Equation 3.11, we get 

{r} = m{a^} 

where 


(3.33) 


m = [L] [iVl 


(3.34) 


B® is the strain displacement matrix in terms of shape functions and is 
given as 


m 


■ dNi_ 
dx 

0 

0 

0 

0 

0 

dy 

0 

0 

0 

dK 

dy 

dN, 

dx 

0 

0 

0 

0 

0 

0 

dN^ 

dx 

0 

0 

0 

0 

0 

Ml 

dy 

0 

0 

0 

Ml 

dy 

dN, 

dx 

0 

0 

dNr 

dx 

N, 

0 

0 

0 

dN^ 

dy 

0 

N, _ 


(3.35) 


Equation 3.30 is divided into a finite number of element domains, Ag, 
within the total domain A in such a way, that 

r E / W\ P\ 

J Ae 

where Ng is the number of elements in the mesh 
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Substituting Equations 3.32 and 3.33 into the above equation, we get 
rt 2 r 

/„ + * = 0 (3.37) 

Le=l 

where 

[M"] = j [iV®]^ [m] [N^] dA 


= U /) M [W'l det [J] dfdt] 

(3.38) 

m = / [B‘f [d] (S'l dA 

•J Ae 


= [B^f [d] [B^] det [J] d^dr] 

(3.39) 

where det[J] is the determinant of the standard Jacobian matrix and 
(■0, T]) are natural co-ordinates of the master element. 

Finally, the element stiffness and mass matrices are assembled to get the 
global stiffness matrix [K] and mass matrix [M] as given below: 

[K\ = £ 1' [bY [n] IB] dxdy 

(3.40) 

and 

[M] ^ [-^1^ [m] [N] dxdy 

(3.41) 


Using Hamilton’s principle, we get the governing equation for free vibra- 
tion of the plate as 


[is:] {a} - [M] {a} = 0 (3.42) 

Assuming harmonic oscillation, 

o = A sin cjt 
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Using the above equation, Equation 3.42 can be written as. 

[K] {a} + [M] {a} = 0 (3.43) 

where {o} is the nodal displacement vector. 

With the above formulation, analysis was done in MSC PATRAN package 
for the present problem. 

The plate was clamped at both the edges along length-wise direction. 
Clamping was done with the help of “C” clamps. Analysis was done to get 
the first three normal modes of the plate and shown in Table 3.1. 

Fig.3.2 shows the discretisation of the plate with QUAD4 isoparametric 
elements. The number of elements are 160 and number of nodes are 189. 
Each node has five degrees of freedom. 

Fig.3.3 to Fig.3.5 shows the first, second and third modes of the plate 
under free vibration. The number 123456 in all these figures at the corners 
indicates that the translations (Tl, T2, T3) and rotations (Rl, R2, R3) are 
constrained along width-wise edges. 

Fig.3.3 shows the first mode shape along length-wise direction. The max- 
imum deformation is at node 11 and the value is 4.79 mm. 

Fig.3.4 shows the second mode shape along the y-direction. The maxi- 
mum deformation is at node 11 and the value is 7.79 mm. 

Fig.3.5 shows the third mode shape along the x-direction. The maximum 
deformation is at node 7 and the value is 4.73 mm. 
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Table 3.1: REAL E IGEN 
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Figure 3.3: FIRST MODE OF PLATE 


40 





Figure 3.4: SECOND MODE OF PLATE 
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Figure 3.5: THIRD MODE OF PLATE 
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Chapter 4 


FORMULATION OF BEAM 
VIBRATION RESPONSE 


The objective of this analysis is to predict the displacement and acceleration 
response of the composite beam embedded with piezoelectric actuator. A 
cantilever beam shown in fig.4.1 was considered for this analysis. Crawley 
and Javier de Luis [ 7] presented one dimensional analysis for a cantilever 
beam embedded with piezoelectric actuators of size equal to beam width 
by considering the strain compatibility in the longitudinal direction. The 
present analysis extends this analysis by enforcing the strain compatibility 
in all directions and also considers finite size of the actuators. 

Certain assumptions which were made for the theoretical formulation are:- 

1. It was assumed that the actuators were perfectly bonded to the sur- 
rounding substructure and the strains in the actuators are compatible 
with the surrounding structure. 

2. A linearly varying strain was assumed across the thickness of beam and 
actuator. 

3. The strains developed in the crystals are within the elastic element. 

Two actuators are placed at equal distances from neutral axis of the beam 
and the direction of the applied electric field along the thickness direction. 
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Sensor 


F 



where 

L = length of beam 
r = location of sensor 
Ij^ = location of actuators 

b = width of beam 
t = thickness of beam 
= width of crystal 

= thickness of crystal 

distance of mid-plane of the crystal 
from neutral axis of the beam 

F = Applied force 

Figure 4.1: CANTILEVER BEAM USED FOR DYNAMIC TEST 
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The beam was assumed to be a linear elastic structure so that the externally 
applied forces could be related to the acceleration, velocity or displacement 
vectors by means of linear influence coefficients. Any continuous structure is 
a system with infinite degrees of freedom, but to avoid the complicacies in 
mathematical modelling, it can be considered as a system with N degrees of 
freedom. 

The lateral deflection v{x,t) of any N degrees of freedom system can be 
approximated by the relation 


v{x, t) = qi{t)(pi{x) + q 2 (t)(p 2 {x) + + qNit)V’N{^) (4-1) 

where g, (z = 1,2, ,N) are generalised co-ordinates and (pi (z = 

1.2, ,N) are assumed dimensionless shape functions which satisfy the 

prescribed geometric boundary conditions. For the cantilever beam, the ge- 
ometric boundary conditions are 


v{x,t) = ^ = 0 (4-2) 

at the built-in end. 

For the cantilever beam, (ft is assumed to equal to 1 — cos((2z — 1) ^). 
Lagrange’s equations of motion are 


^\dqi) ^q^ dqi 


(4.3) 


where T =Total kinetic energy 
V = Total potential energy 

= 1, 2, ...., iV) = generalised forcing functions 
— q\iQ2i ) generalised coordinates respectively. 


Thus 





(4.4) 


The total kinetic energy of the beam will be sum of the kinetic energy 
of the beam and kinetic energy of the sensor or accelerometer located at a 
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distance I from the built-in end. 
Kinetic energy of the beam 


Ti 



(4.5) 


where m — mass per unit length of the beam. Kinetic energy of the sensor 
is 


at X = / , where rris is the mass of the sensor. 

Total kinetic energy 

T = Ti+T2 (4.7) 

Substituting Eq.4.4 into Eq.4.5 and Eq.4.6, 



(4.8) 


^2 = hris fl - cos ttz) fl - cos ) ttI ) (4.9) 


^ T=i j=i '• ^ 2L 

Evaluating the above terms, we get 
^ ^ .m r, 2L 


^ ^ m ( 
1=1 7=1 I 


2L 


or r 


m, 

"2 


( 


2z- 1 
2L 


* ' ^ 1 - cos f ------ ) ttZ ) ( 1 - cos ( ttI^ |] (4.10) 


2i- r 

2L 


The total potential energy V will be 



dx 
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N N 



Trx cos 



TTXdx 


v = 


EbIb'^^ 

64L3 


N 




t=l 


(4.11) 


To employ the Lagrange s principle, we need to know the virtual work 

done by the externally applied force T’ at the tip of the beam. This can be 
calculated as 






tvL 


N 


5W = F5qi 
1=1 


2i-V 

2L 


)’rsin (^^)7rl, (4.12) 


the above expression of 5W can be written as 


N 


sw = '£QM, 

Z=1 


where 


Qi = Fi 


2i-l 


TTsm 


2z- 1 


ttL 


V 2L j ~-\ 2L 
Applying the Lagrange’s equation to obtain the equation’s of motion 


(4.13) 


^ _ 

dt ^q^ 
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The term is deleted as there is no term in kinetic energy equation 
which is a function of q^. 


m( 2L 2L 

( ^ - 2_^ L 1 -C- + 77: — TT“(-1)* + 7 + 5. 

^=lJ=l 2 [ {2i-l)7r^ ’ (2j_l)7r^ ^ ^ ‘ 


L' 

2 


y r- ( 


2 


{( 


22-1 




^ ^ ^ ( 03 

1=1 j=i L (2? — 


l)7r< + 




(4.14) 


dfdT\^^mj 2L , ,, , 2L L] 


m. 


{ i%r) {'■ - (^) ”■') }i+ 


tt I? " ( 2 ^<-^)‘ - - ^ 4 ) + 


4 { (1 - {^) 4 (1 - {% r ) 4 }i 

The above term in the right hand side can be written as My$i where 
are the elements of the mass matrix whose values are given by 


M 


y 


3 4 , .,il , /22-r '2- 

m{- - — :t— (-1) [> + rus j 1 - cos 


2 (22 - l)7r 


2L 


-) TT^y 


^fi=j 

(4.16) 
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2 


M, 




.m 


1 + 


(2z - 1) 


TT 


{~ir+ 


{2j - 1) TT 


(-IF + 


The mass matrix is a symmetric square matrix of the order n x n and the 
off-diagonal terms are non zero, is a vector of n x 1 order which can be 
written as 


= {91,92, ,9n} 


(4.18) 


and 


[M] = 


Mu 

M21 M22 


Afnl 



(4.19) 


Similarily, the potential energy term differentiated with respect to q in 
the Lagrange’s equation will give rise the stiffness matrix in which elements 
are given by the term 


K., = {2i - l)U.i (4.20) 

where 6ij is the Kronecker delta function. Substitution in terms of matrix 
will give the equation as 


m {?} + m { 9 } = { 0 } ( 4 . 21 ) 

Natural frequencies of the beam can be calculated from the above equa- 
tion by equating the forcing function to zero. So, for evaluating the natural 
frequencies, the equation will be 
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Substituting 


M {9 } + m {9} = { 0 } 


(4.22) 


{q} = {A} sin {ujot + (4.23) 

in above equation, gives 



where {0} is a n x 1 order zero vector. 

After premultiplying the above equation by [K]~^ , it can be written as 



We now define a column vector 



which is the eigen values vector of the matrix [AT] 

Corresponding to each of the eigen values, there will be an eigen vector. 
In order to decouple the equations of motion, the eigen vectors must be 
normalized. The normalizing procedure, which has been used here involves 
adjusting each modal amplitude which satisfies the condition 
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{Z,f [Jf] {Zr] = 1 (4.24) 

where Z^ = ■- {/I)} and (A)} is the eigen vector corresponding to the 
eigen value ^ 

Corresponding to each eigen value of the set = 1,2,3, ,n) 

there will be a normalized eigen vector = 1, 2, 3, ,n). 

cl can be calculated for each eigen vector by making the substitution 

A consequence of this type of normalizing together with the modal or- 
thogonality relationships relative to the stiffness matrix is 

= I (4.25) 

where $ is the complete set of N normalized mode shape vectors and I 
is a n X n identity matrix. The mode shapes normalized in this way are said 
to be orthonormal relative to the stiffness matrix. 

Using the normalized vector, the equation 

m {A} = [iq {A} (i) 

can be written as 

lM]{Zi} = lK][^^ 

similar equations hold good for w<, 2 ,Wo 3 , with corresponding 

eigen vectors {A^} , {Af} , , {A”} . In the matrix form, it can be writ- 

ten as 

[M] [zl z? z;] = [AT] 

or 

[M] [ZJ = [K] [^] m 


^ol ^02 




U)\ 


on J 


(4.26) 
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«rncfK v.U-rf 

g?o a... 1.3.7 9.Q.^ 



where 


[^o] 


ol 


Premultiplying Eq.4.26 by [Zf gives 

[Zf [M][Z] = [zf [K][Z][^o] 

Substituting Eq.4.24 in above equation gives 

= (4.27) 

Introducing the transformation {q} — [Z] {p} in the Eq.4.21 gives 


m [z] {9} + [K] [z] w = { 0 } 


(4.28) 


where {Q} — < 


Qi 

Qi 

Qn 


is forcing function vector 


we have Qt = F sin and forcing function F = Fg sin ut 

Premultiplying Eq.4.28 by [Z]^ 


[Zf [M] IZ] {9} + [Zf [K] IZ] {p} = [Zf {Q} 

Substituting Eq.4.24 and Eq.4.27 in the above equation, reduces to 

I'i'J {9} + {?} = [Z'f {«} (4-29) 

The equations are decoupled, the solution is 
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1 


p^ = 



Zjt {Q} smut 


j=i 


(4.30) 


Substituting back for {9} — [Z] {p} , we can calculate the displacement 
and acceleration in terms of generalised coordinate. 


Experimental Procedure For Dynamic Test 

Manufacture of glass/epoxy specimen 

The smart glass/epoxy test specimen was prepared using plain woven glass 
fabric material. The dimensions of the test specimen, actuators and sensors 
are shown in Fig.4.1 and Table 4.1. 

Embedding of PZT crystals 

Leads of diameter 0.2 mm were flattened and attached to crystals on both 
faces by soldering. The overall thickness of the crystal along with leads is 
about 0.40 mm. So, to accommodate these crystals in the laminate, fibers 
from three lamina’s were removed using a blade. Holes were 1 mm larger than 
the dimension of the PZT crystals and the crystals were inserted without 
causing any change in the nominal thickness of the laminate. Finally, the 
laminate was first cured at room temperature for 12 - 18 hours and then 
hot curing was followed. Aluminium end tabs were fixed to the beam at the 
clamping end on top and bottom to avoid delamination due to clamping. 

Experimental Procedure for Dynamic Test 

The specimen was clamped horizontally between steel plates. The base plate 
was welded to the fixed structure, specimen was placed on it and the top 
plate was tightened over the beam by bolts, which was effective in providing 
a cantilever boundary conditions. The natural frequency of the laminate 
was measured with the help of FFT Analyser for the first three modes and 
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Table 4. 

L: DYNAMIC TEST SPECIMEN DETAILS 


Material 

Dimensions (mm) 

Laminate 

glass/epoxy 

260 X 80 x 1.61 

Crystal 

PZT 

10 X 10 X 0.25 

Sensor 

Endevco 2242C Accelerometer 

Sensitivity 2.22 pc/g 


were compared with the analytically predicted value. Analytically predicted 
values are shown in Table 4.2. 

The beam is connected to a magnetic shaker through a connecting rod 
at the tip of the beam. The beam is subjected to concentrated time varying 
load. The forcing frequency was 65 Hz. The excitations of the beam were 
sensed by the embedded crystals. The output of the crystals was amplified 
by using a charge amplifier and measured with the help of FFT Analyser. 
The experimental set up is shown in Fig.4.2. The response of the PZT crystal 
and sensor is shown in Figs.4.3 &: 4.4. 
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Table 4.2; REAL EIGEN VATJJES OF TH E CANTILEVER BEAM 
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Where 

A; PZT Crystals 
B : Accelerometer 

C : Charge Amplifier 

D : FFT Analyser 

E : Magnetic Shaker 

Figure 4.2: SCHEMATIC DIAGRAM OF EXPERIMENTAL SET UP 
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Figure 4.3: FORCED VIBRATION RESPONSE OF BEAM BY ACTUATOR 
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Figure 4 4: FORCED VIBRATION RESPONSE 


OF BEAM BY SENSOR 
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MSC/PATRAN Version 9.0 03-Fel3-02 13:00:11 

Fringe: DefaoH, Mode 1:Freq.=12.199: Eigenvectors, Translational-(NON-LAYERED)(ZZ) 
Deform: Default, Mode 1:Freq.=12.199: Eigenvectors, Translational (NON LAYERED) 
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MSC/PATRAN Version 9.0 03-Feb-02 13:01 :55 

Fringe; DefaoH, Mode 2;Freq.=76.433: Eigenvectors, Translationa!-(NON-LAYERED) (ZZ) 8 52+00 

Deform: Default, Mode 2:Freq.»76.433: Eigenvectors, Translational 

7 55+00 



Figure 4.7: SECOND MODE OF BEAM 
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MSC/PATRAN Version 9 0 03-Feb-02 13:03:15 

Fringe: Defaali, Mode 3:Froq.-2 1 3.93; Eigenvectors, Transiational (NON-lJ^yERED) (ZZ) 5 61 +00 

Deform: Defaatt. Mode 3:Froq.»213.93; Eigenvectors, Translational 

4 67+00 



Figure 4.8: THIRD MODE OF BEAM 
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Chapter 5 

RESULTS AND CONCLUSIONS 


5.1 Composite Laminate 

The composite laminates were prepared as described earlier. They were 
characterized for mechanical properties. The tensile test was carried out on 
a 10 ton MTS servo hydraulic machine. The Elastic modulus was found to 
be 5.776 GPa experimentally.The Burn out test was carried out to calculate 
the exact fiber volume fraction. The fiber volume fraction was found to be 
0.46 and the rule of mixture was applied to calculate the Elastic modulus 
of the laminate. The value obtained by this method was 35.01 GPa. The 
poisson’s ratio was found to be 0.22. Tensile tests were performed for two 
specimens on a table top tensile testing machine. The modulus was found to 
be 5.50 GPa. 

The procedure suggested by Naik[22] was adopted to determine the mod- 
uli. The lamina was observed through a optical microscope to get the yarn 
width and gap width. Based on this information, the fiber volume fraction 
and poiason’s ratio, the Elastic modulus and Shear modulus was found to be 
16.50 GPa and 3.60 GPa respectively. 
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5.2 Dynamic Response 

Natural frequency of the clamped - clamped plate and cantilever beam was 
measured theoretically and experimentally. Both the values for first three 
modes of the plate and beam are given in Tables 5.1 and 5.2 respectively. 
The experimental and theoretical acceleration of the sensor is shown in Table 
5.3. The theoretical frequency modes of the plate and beam are shown in 
Figs. 3.3 to 4.8. 

5.3 Conclusions 

MSC Patran was used for the analysis of both the plate and beam problem. 
The plate was meshed using Quad4 elements as shown in Fig.3.2. The anal- 
ysis was carried for normal modes of the plate. The beam was discretised 
using Bar2 elements as shown in Fig.4.5. Sinusoidal loading {Fosinut) was 
applied at the tip of the beam. The forcing frequency was 65 Hz. Analysis 
was done and the theoretical acceleration was found. 


5.4 Future Work 

With the experience gained in this work, one can address more complex 
structures such as plates and shells and can optimize for placements of ac- 
tuators and sensors. The studies can be extended to close loop behaviour, 
where the response of the structure can be controlled through the complete 
cycle of sensing, processing of information and actuation. 
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Table 5.1: COMPARISON OF FIXED-FIXED PLATE NATURAL FRE- 
QUENCIES 


Mode 

Experimentally Predicted 
Frequency (Hz) 

Theoretically Predicted 
Frequency (Hz) 



166.36 

Second 

325.00 

310.64 

Third 

475 00 

461.41 


Table 5.2: COMPARISON OF CANTILEVER BEAM NATURAL FRE- 


Mode 

Experimentally Predicted 
Frequency (Hz) 

Theoretically Predicted 
Frequency (Hz) 

First 

22.50 1 

12.20 

Second 

135.00 

76.43 

Third 

275.00 

213.93 


Table 5.3: COMPARl 

[SON OF SENSOR ACCEL] 

ATTON IN THE BEAM 

Applied Force (N) 

Experimentally Predicted 

Theoretically Predicted 


Acceleration (m/sec^) 

Acceleration (m/sec^) 

50 

52.024 

65.250 
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